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ABSTRACT 

We investigate the contributions coming from the penguin operators in 
the nonresonant B~ —>■ MMtt~ (M = ir~,K~,K°) decays. The effective 
Wilson coefficients of the the strong penguin operators 04,6 are found to be 
relatively larger. We calculate the contributions arising from the O4 and 0$ 
operators in the nonresonant decays B~ — > MMn~ (M = n~, K~, K°) using 
a model combining heavy quark symmetry and the chiral symmetry, devel- 
oped previously. We find that the forbidden nonresonant B~ — > K K°tt~ 
decay occurs through the strong penguin operators. These penguin contri- 
butions affect the branching ratios for B~ — > M M tt~ (M = 7r~, K~) by only 
a few percent. The branching ratio for B~ — > K°K°tt~ is estimated to be of 
the order 10 -6 . 



There is considerable interest in understanding the decay mechanism of 
the nonleptonic charmless three body decays of B mesons |3]> @> El- The 
importance of penguin operators in three body decays of charged B mesons 
has recently been questioned Jl|. In the analysis of the Dalitz plot for B~ — > 
7r + 7r~7r~ the authors of [Q have assumed that the nonresonant decay ampli- 
tude is flat, having no dependence on the Dalitz variables. They also assumed 
that the contributions of the penguin operators can amount to as much as 
20% of the dominant decay amplitude. Others have made predictions for the 
branching ratios of decays []2], [3], ||] motivated in part by the CLEO limits on 
some of the nonresonant decays of the type B + — > h + h + h~ |J. CLEO found 
the upper limits on the branching ratios BK(B + — > 7t + 7t~tt + ) < 4.1 x 10~ 5 
and BR(B + -> K+K--K+) < 7.5 x 10 -5 . In addition there is hope that the 
CP violating phase 7 can be measured from the asymmetry in charged B 
meson charmless three body decays |], |5j . 

Motivated by the need to understand whether the nonresonant decay 
amplitudes for B~ — > MMn~ (M = tt + ,K + ) involve significant effects due 
to the penguin operators we have investigated the contributions coming from 
the penguin operators - |12] in these nonresonant decay amplitudes. The 
decay B~ — > K°K°it~ is CKM forbidden [13, 14]. However, we found that 
B~ — > K K 7T~ occurs through penguin operators. A measurements of this 
rate would allow one to extract the product of the CKM matrix elements 

v tb v; d . 

In our analysis we will use of the factorization approximation in which the 
main contribution to the nonresonant B~ — > M Mtt~ amplitudes come from 
either the product < MM\(ub) v ^ A \B~ > < ir-\(du) v _ A \0 > 
or < 7r~ M\(db)v-A\B~ > < M\(uu)v-a\0 >, where (qiq2)v-A denotes 
gi7 M (l— 7 5 )g 2 - For the calculation of the matrix element < MM\(ub)v-A\B~ > 
we extend the results obtained in |Tj|, where the nonresonant D + — > K~tt + Iv 
decay was analyzed. The experimental result for the branching ratio of the 
nonresonant D + K~n + lv decay was successfully reproduced within a hy- 
brid framework [JT5] which combines the heavy quark effective theory (HQET) 
and the chiral Lagrangian (CHPT) approach. The combination of heavy 
quark symmetry and chiral symmetry has also been quite successful in other 
analyses of D meson semileptonic decays |16 1 - [^] . 

Heavy quark symmetry is expected to be even better for the heavier B 
mesons |L9|, ^0[ . However, CHPT might be less reliable in B decays due to the 
large energies of light mesons in the final state. It is really only known that 
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the combination of HQET and CHPT is valid at small recoil momentum. 
To take into account the larger recoil energies of the light mesons in our 
previous work PH ], we modified the hybrid model of |TB[ - [^tj to describe 



the semileptonic decays of D mesons into one light vector or pseudoscalar 
meson. Our modification is quite straightforward: we retain the usual HQET 
Feynman rules for the vertices near and outside the zero-recoil region, but we 
include the complete propagators instead of using the usual HQET propagator. 
This quite reasonable modification of the hybrid HQET and CHPT model 
enabled us to use it successfully over the entire kinematic region of the D 
meson weak decays |L5], ^1], [22]j . 

In the following we systematically use this model to calculate the contri- 
butions of the penguin operators to the nonresonant B~ — > MMir~ (M = 
it~,K~,K°) decay amplitudes. We first analyze the contributions coming 
from the 46 penguin operators [12], since their effective Wilson coeffi- 
cients are the largest. We then determine the dependence on the Dalitz plot 
variables. The operator O4, as defined in P, [12| , has the same dependence 
on the Dalitz plot variables as the tree- level operator 0±, while Oq exhibits 
different energy dependence. Finally, we discuss the influence of these oper- 
ators on the branching ratios for B~ — > M M~ (M = 7r + , K + ) and estimate 
the branching rate for B~ — > K K°tt~. 

The effective weak Hamiltonian for the nonleptonic Cabibbo-suppressed 
B meson decays is given by [§, 0, |T2| 



H eff = ^[v: d v ub (c 1 o lu + c 2 o 2u ) + v; d v cb ( Cl o lc + c 2 o 2c )} 

10 

- ^([K&KX + v cbV* d c1 + V tb V t a4)Oi] + h.c. (1) 

i=3 

where the superscripts u, c, t denote the internal quark. The operators Oj are 
defined in [§, [K], 0. We rewrite 3 — 6 , using the Fierz transformations, 
as follows: 

3 = E d lfl (l- l5 )bqY(l-l,)q, (2) 

q=u,d,s,c,b 

Oa= E d7M(l-75)gg7 M (l-75)&, (3) 

q=u,d,s,c,b 

5 = E d^(l- l5 )bqY(l + l 5 )q, (4) 

q=u,d,s,c,b 
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6 = -2 £ d(l- l5 )qq(l + l5 )b, (5) 

We consider two possibilities for the Wolfenstein parameters p and r\: Case I, 
the effective Wilson coefficients a\ * determined in M for b — > d transitions 
with iV c = 3 , and p = 0.12 and rj = 0.34: 

af f = -0.0412 - 0.0036i (6) 

and 

a e / f = -0.0548 - 0.0036i. (7) 

Case II, the effective Wilson coefficients a e / determined in [|§ for b — > d 
transitions with iV c = 3 and p = 0.05 and r\ = 0.45: 

af * = -0.048 - 0.007z (8) 

and 

a e 6 ff = -0.060 - 0.007i (9) 

The effective coefficients al and a§ are one order of magnitude smaller 
[§, U than these two and therefore we can safely neglect the contributions 
arising from 3 5 operators. 

The quark currents required in the weak Hamiltonian ([]]) can be expressed 
in terms of the meson fields, as previously described explicitly in ||, 13, 



The operator Oq can be rewritten as the product of the density operators. 
For the d(l — 75)? scalar and pseudoscalar quark density operator we use the 



CHPT result [23]. The explicit chiral symmetry breaking, to lowest order in 



the chiral expansion, is obtained by adding the quark mass term p3 



f 2 

£ s = ^-{trB(MU^ + UM^)}, (10) 

where M = diag(m u , m^, m s ) and B is a real constant that can be expressed 
in terms of quark and meson masses; e.g., to lowest order m 2 R0 = B(m s + m^) 
and U = exp(i2U/ f) where II is a pseudoscalar meson matrix. Using (|lTi|) 
one can easily bosonize the density operators: 

ft(l - 75)^ = BU 3i (11) 
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For the calculation of the density operator q{l + 75)6 we use the relations 



q l5 b = — a (g7 a 75&), 
m b 



and 



qb 



m b 



daiq^b), 



(12) 



(13) 



where m q has been dropped since m q « m b . 

The evaluation of the matrix elements < M\q(l + 7 5 )6|5 > and < 
MM\q(l + 7 5 )&|i? > can then be reduced to the evaluation of the matrix 
elements of the weak currents < MM | (77^756 |i? > and < M\q^^b\B >. As- 
suming factorization, we evaluate the matrix elements of the operator 0§ : 

<MMM\0 6 \B> = -2 ]T {< M\d(l - 7 5 )g|0 >< MM\q{l + ^ 5 )b\B > 

u,d,s,c,b 

+ < MM| J(l - 75 )g|0 >< M\q(l+j 5 )b\B > 

+ < MMM\d(l--f 5 )q\0><0\q(l + -f 5 )b\B >}. (14) 

The matrix elements < M\d(l - 75)g|0 >, < MM\d(l - 7s)g|0 >, and 
< MMM\d(l — 75)g|0 > are easily calculated using (|i~l"|). For the calculation 
of the matrix elements < M]g7 M &|S > and < MM|g7yy5&|-6 > we generalize 
the results obtained in the analysis of D meson semileptonic decays described 
in detail in [15] and [21]. The matrix element < M|g7 jU (l — 75)6]!? > is given 
by |1], H 



<M(p')\q lf ,(l- l5 )b\B(p B )> = [(pb+p). 



in 



H 



in 



M 



+ 



rn 



H 



rn 



q»F {q 2 



(15) 



where q = Pb — p' and -Fi(O) = -Fo(O). The form factors are found to be 
m. I22H 



and 



2 V m B'* 



+ — 2~ 7T+9fB'*\ , (17) 

where B'* denotes the relevant vector meson pole and g is the B*BM coupling 
constant. 

To evaluate the matrix element < Mi(pi)M 2 (p 2 )\(qib)v-A\B~ (p B ) > we 
will also use and generalize the results obtained previously in the analysis 
of the nonresonant D + — > -k + K~Ivi decay width Hl5| . We write the matrix 
element < Mi{pi)M 2 {p 2 )\{ub)v~A\B~ {p B ) > in the general form 

< Mi(pi)M 2 (p 2 )|«7^(l -ls)b\B~{p B ) > = ir{p B -p 2 -p 1 ) fl 

+iw + (p 2 + px)^ + iw^(p 2 - pi)^ - 2he m ^p a B p 2 p\ . (18) 



The form factors wj r for the nonresonant decay are given in 0, [15 



<(pi,p 2 ) = -JL S * m % m *\ {l Pi-(p B - P i) 



fih (pb - Pi) 2 - m Bif 2m 2 m 
h __ Vrn B a 2 1 
2/i/ 2 2/ x / 2 m| J 



JB y/HL B LX 2 L 

+ TTT - ? n * 2 Pb • (P2 - Pi), (19) 



W-(Pl,P2) 



+ 



3/2 1/2 

/1/2 (pb 



Pi) 



■[1 + 



Pi • (P 



B 



Pi 



B* 



2mL 



(20) 



The parameters «i j2 are defined in Note that both the a± and a 2 terms 
are important in fll9"|) and (p0|), which was previously overlooked 0. Within 
this same framework [15|, |2l| we evaluate r nr 



\Pl,P2) 



l+{3 



hh 



Pb ■ (P2-P1) 



\m B > 



h 



1% H ffB"m* B ,m B t 



m B (p B - pi 
1 



P2j 



— m 



B' 



m B 



hh 



(Pb ~ Pi - P2) 2 - my, 



5 



X 



+ 



[Pi ■ P2 ~ t^P2 ■ (Pb - Pi)pi ■ (pb - Pi)} 

(Pb - Pi) 2 - m 2 Bljf 
2g f B >*m B ,* Pi • (p B - Pi) 
hh (Pb-Pi) 2 - m B ,, m| v 
f B . ot 2 \[m 



+ ^f T + ^^i R PB-{P2-Pi) (21) 



2/1/2 2m 



Here B', B *, B" denote the relevant B meson poles, and f\ y2 denotes the 
pseudoscalar meson decay constants. The coupling (3 has been analyzed in 
fm and found to be close to zero and therefore will be neglected. 

The matrix element of the operator O4 can be evaluated straightforwardly 
using factorization: 

< 7r _ 7r + 7r _ |0 4 | J B _ > nr = < ii + ix~ |«7 M (1 - 1<o)b\B~ > nr 

<7r-|J 7 ^(l-75)n|0>, (22) 

and the corresponding expression for the B~ — > K~K + 7t~ matrix elements 
can simply be obtained by the replacement ir + and tt~ by K + and K~ respec- 
tively. Note that the matrix element < MM\q% 7 M (1 — 75)9210 > is dominated 
by resonant contributions. Using the variables s = (p B — P3) 2 = [p2 +P1) 2 , 
t — (.Pb — Pi) 2 = G°2 +P3) 2 u = (pb — P2) 2 = (pi + P3) 2 and the pseudoscalar 
meson decay constants fa, we can then write the nonresonant decay matrix 
element of 4 as 

<M 1 ( Pl )M 2 (p 2 )n-(p 3 )\0,\B-(p B )> nr = [f„mlr nr (s,t) + ^(mi 

-s-ml)w n + r (s,t) + ^(2t + s-m 2 B -ml-m 2 2 -ml)w m Xt)} 

(23) 

where Mi, M 2 represent either 7r + , 7r _ or K + , K~ . 

Using factorization the matrix elements of 0§ can be written as 

< ■K + (px)'K~(P2)'K~{Pz)\0 & \B~(pb) >„.r = 

-2{< n + {p 1 )7i-(p 2 )\u-f fl b\B-(p B ) > X 
n 

<7r-(p 3 )|rf(l-7 5 )w|0> + 

m B 
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< n + (pi)n (p 2 )n (p 3 )|d(l-75)«|0> X 

<0\u^b\B-(p B )>^ + 
m B 

< 7i + ( Pl )7i-(p 2 )\d(l - 7 5 )w|0 >< 7r~(p 3 ) \u^b\B-(p B ) > x 

P -^^ + ( Pl ~ p 3 )h (24) 
m B 

where we have assumed m;, ~ m B . The corresponding result for B~ — > 
7i~K + K~ can be straightforwardly obtained simply by replacing tc + tt~ by 
Using the expressions for the matrix elements of the current and the 
density operators we find 

B 

<M 1 (p 1 )M 2 (p 2 )7r-(p 3 )\0 6 \B(p B )> nr = -/ 3 

m B 

{-[m\r nr + (m B - m\ - s) + w r "\2t + s - m B - m\ - m\ - ml)w™ r } 

+[(m| - ml)F Q (s)] - \^m B ). (25) 

For the B~ — > 7r~7r + 7r~ decay there is an additional term with the replace- 
ment s t, since there are two identical pions in the final state in this 
case. 

The nonresonant amplitudes for the B~ — > MMtt~ (M = 7r~, K~) decays 
can be written in terms of the following matrix elements 

M nr (B- - MMn-) = ^{V ub V: d a\ ff KMM-u-p^B- > 

+V tb V t * d (a/ f < MMtt-\0 4 \B- > + al n < MMtt-\O q \B- >)}. (26) 

The matrix element < MMn~\0\\B~ > (M = n~,K~) was analyzed in [|]]. 

Contrary to the CKM-allowed cases in which the main contribution comes 
from the operator Oi, we notice that the CKM- forbidden decay B~ — > 
K K 7t^ occurs through the contributions of the penguin operators O4 and 
Oq. The nonresonant matrix elements are 



<K\p l )K-(p l )K\p 2 )^(pz)\O i \B~{p B )> m . = l^[ m lr nr (s,t) + (m 2 



B 



-t-m 2 K )wl r (s,t) + (2t + s-m 2 B -2m 2 K - m>! r (s,i)] (27) 
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and 



< K\ Pl )K ( P2 )ir-(p 3 )\O 6 \B-( PB ) > nr 
w™'{s, t) [m\-t- m 2 K ) + wl r {s, t){2t + s 

+ [(m 2 B -ml)F (s)} 



B—{r nr (s,t)m 
m B 

ml) 



m\ — 2m K 



4 Uf. 



B 



3 / 



m B } 



(28) 



K 



The nonresonant amplitude for the B (p B ) — > K (pi)K (p 2 )Tc decay can 



be written in terms of these matrix elements (]27|) — (|28|) 

G p . 



M nr (B- -> K K 
{af f < K°K%-\0 A \B- > 



V2 



v tb v; d 



a, 



eff 



< K°K°7T-\O e \B- >}. (29) 



The partial widths for the nonresonant decay B 
ir~ , K~ , K°) is given by 



MM 



7T 



(M 



T nr (B~ 



MMtt- 



(2tt) 3 32m| 



y |-M nr | 2 <is (it. 



(30) 



In the numerical calculation of the branching ratios we follow the discus- 
sion of the input parameters given in [|]]. From heavy quark symmetry we 
have used j B j Jd = \/rn D /m B with the reasonable choice /o ~ 200 MeV 
The B decay constant is then f B ~ 128 MeV. In we found that 



22, 26 



the parameters af p = -0.13 GeV 1 / 2 and a 2 p = -0.36 GeV 1 / 2 lead to the 
branching ratio BR(B - — ► 7r~7r + 7r + ) being smaller than the experimental 
upper limit || and we rejected this possibility. Here we also use the values 
of ai t 2 as in ||]. And, as discussed in [|J, here we also consider the range 
0.2 < ' g < 0.23. 

In our numerical calculations we considered both cases I and II: the effec- 
tive coefficients a e JJ given in @, (^) from Q (case I), and a e JJ given in (|8|), 
(|Sp from (case II). Then for the CKM matrix we must use the correspond- 
ing input parameters in the Wolfenstein parametrization of the CKM matrix 
( V ub = AX 3 (p - ir,), V ud =l- A 2 /2, V td = AX 3 (1 -p- irj), V tb = 1). The 
numerical value of B can be determined from B = (2m 2 K — m 2 )/2m s . Taking 
m s (p = 5 GeV) = 150 MeV, the same value used in || for the extraction 
of the effective Wilson coefficients, we find B — 1.6 GeV. Inspection of the 
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contributions coming from the 04,6 operators shows some cancellations occur 
among the combinations a% O4 and a% 0§. One also can explicitly see the 
dependence on the Dalitz variables. 

In Table I we present the penguin contributions of the operators to 
branching ratios for the B~ — > 7r~7r + 7r - and B~ — > K~K + tt~ together with 
the dominant tree level contribution of the operator 0\. Both numerical 
results I and II as well as the range of g, as discussed above are presented. 

It is clearly evident from these quantitative numerical results that the 
uncertainties coming from the input parameters give much larger uncertain- 
ties in the branching ratios than the contributions of the penguin opera- 
tors. Interestingly the penguin contributions, while small, are less sensitive 
to the input parameters than the dominant tree level contributions, which 
is quite sensitive to these input parameters. Since the amplitudes for the 
B~ — ► 7r + 7r~7r + and B~ — > K + K~tt + decays receive rather small corrections 
from the penguin operators we do not expect significant changes in the CP 
violating asymmetry, which we have discussed in ||. 

We also calculated the branching ratio for the forbidden nonresonant 
decay B~ — > K K°ir~ finding 

8.4 x KT 7 < BR(B~ -> K K\-) < 8.7 x KT 7 , (31) 

for the Case I (see @ and (0)) and the range 0.2 < g < 0.23. For the Case 
II (see (8) and (9)) we found 

1.4 x 10~ 6 < BR(B~ -> K°K°tt-) < 1.5 x 10" 6 , (32) 

for the range 0.2 < g < 0.23. Measurement of this branching ratio is impor- 
tant as it provides information about the effective Wilson coefficients all . ^ 
is interesting to note that in the factorization approximation, as mentioned 
earlier, this decay is entirely induced by the penguin interactions. Final state 
interactions (FSI) effects, could alter this, however we believe this is unlikely 
as data on color-suppressed decays indicate that the branching ratio for B° 
to D° and a neutral light hadron is indeed suppressed. Thus we expect FSI 
would contribute at most a branching ratio for K°K° mode of the same order 
as the penguin terms. This could be checked in future measurements of this 
decay rate. 

To summarize, we have quantitatively analyzed the penguin contributions 
to the nonresonant B~ — > MMir~ decay amplitudes (M = 7r~, K~ , K°), in- 
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eluding the dependence on the Dalitz variables. We calculated the branching 
ratios for B~ — > MMn" decays (M = n~,K~) including the penguin con- 
tributions and found that they can possibly change the branching ratio as 
much as 15%. However, while the penguin contributions are small and not 
very sensitive to the uncertainties in the input parameters, the corresponding 
uncertainties in the dominant tree level contributions are considerably larger 
than the penguin contributions. We also found that the branching ratio for 
the CKM forbidden nonresonant decay B~ — > K K°ix~ is of the order 10~ 6 
and is entirely induced by penguin effects. 

This work was supported in part by the Ministry of Science and Tech- 
nology of the Republic of Slovenia (S.F), and by the U.S. Department of 
Energy, Division of High Energy Physics under grant No. DE-FG02-91- 
ER4086 (R.J.O.). S.F. thanks the Department of Physics and Astronomy 
at Northwestern University for warm hospitality during her stay there where 
part of this work has been done. 
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p 


V 


9 


BR(B~ -> n-n+ir-) 


BR(B~ -> K+K-7V-) 


0.12 


0.34 


0.2 


(3.1-0.4) x 10~ 5 


(5.3-0.1) x 10~ 5 


0.12 


0.34 


0.23 


(3.5-0.5) x 10~ 5 


(5.8.-0.1) x 10~ 5 


0.05 


0.45 


0.2 


(4.9 - 0.4) x 10~ 5 


(8.7-0.1) x 10~ 5 


0.05 


0.45 


0.23 


(5.5-0.4) x 10~ 5 


(9.0-0.03) x 10~ 5 



Table 1: The branching ratios for B~ -> MMtc~, (M = -k'^K') for 
two cases of p and rj parameters determined in [8] and [9], described in the 
text, respectively. The first number in the brackets is the main contribution 
coming from the operator 0± and the second number is the contribution to 
the branching ratio coming from the operators O^q. 
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